Synthesis, crystal structure and magnetism of /3-Fei oo(2)Sei.oo(3) single crystals 
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Understanding iron based superconductors requires high quahty impurity free single crystals. 
So far they have been elusive for /3-FeSe and extraction of intrinsic materials properties has been 
compromised by several magnetic impurity phases. Herein we report synchrotron - clean /3-FeSe 
superconducting single crystals grown via LiCl/CsCl flux method. Phase purity yields evidence for 
a defect induced weak ferromagnetism that coexists with superconductivity below Tc. In contrast 
to Fei+yTe - based superconductors, our results reveal that the interstitial Fe(2) site is not occupied 
and that all contribution to density of states at the Fermi level must come from in-plane Fe(l). 

PACS numbers: 74.70Xa, 61.72y, 74.62Dh, 74.62Bf 



I. INTRODUCTION 

The physics of complex superconductors, such as the 
cuprates and iron-based superconductors, cannot be un- 
derstood unless pure, high quality materials are available 
that allow the intrinsic properties to be separated from 
extrinsic and impurity effects. In FeAs and Fe(Se)Te, 
just as in high-Tc cuprate and heavy fermion materials, 
competing or coexisting magnetic order is closely asso- 
ciated with superconductivity.^ This suggests proximity 
to a magnetic critical point and an unconventional origin 
of superconductivity where spin fluctuations may con- 
tribute to pairing. ^^ The observation of weakly local- 
ized rather than itinerant magnetism sensitive to struc- 
tural changes raises the fundamental question of how 
strongly correlated are the charges in Fe superconduc- 
tors and what is the origin of the magnetic order?^^ 
Of particular interest is superconducting /3-FeSe, a com- 
pensated semimetal, without a crystallographic charge 
reservoir, that super conducts at about 8 K without any 
carrier doping.^ It has a giant pressure coefficient of Tc of 
9.1 K/GPa enhancing Tc up to a maximum of 37 K, the 
third highest known critical temperature for any binary 
compound^ 

A major obstacle in understanding intrinsic magnetism 
in /3-FeSe has been the purity of the material itself. Mag- 
netic impurities such as a-FeSe, FeySeg, Fe304 and ele- 
mental Fe are ubiquitous in all as-grown crystals and 
sometimes poly cry stalsJ^"— They contribute to the large 
ferromagnetic (FM) background, seen in the M-H loops 
below superconducting Tc. Modification of the original 
Fe-Se phase diagram near 1:1 stoichiometry suggested 
that /3-FeSe is not stable at the room temperature since 
it converts to hexagonal a-FeSe below 300° C.^^^^^ Con- 
sequently, the absence of an exposed liquidus surface in 
the binary alloy phase diagram and the metastable na- 



ture of the superconducting FeSe are considered to be 
prohibitive and insurmountable factors for single crystal 
preparation using standard synthesis methods. 

Here we describe a synthetic approach that yields stoi- 
chiometric and phase pure material and we report intrin- 
sic structural and magnetic properties of superconduct- 
ing /3-FeSe. These include evidence for defect - induced 
weak ferromagnetism (WFM) and the absence of inter- 
stitial Fe(2) whose occupancy governs the magnetic and 
structural phase diagram in isostructural Fei+^Te^ 



II. EXPERIMENT 

Powders of LiCl and CsCl, elemental Fe and Se were 
added together with the flux into an alumina crucible and 
sealed under partial Ar atmosphere. The ampoule was 
heated to a homogenization temperature of 715°C, where 
it was kept for 1 h and then removed into a preheated 
furnace at 457° C. After slow cooling to 300° C, it was 
quenched in water. 

Medium resolution, room temperature (300 K) X-ray 
diffraction measurements were carried out at X7B beam- 
line at National Synchrotron Light Source (NSLS) at 
the Brookhaven National Laboratory, using a 0.5 mm^ 
monochromatic beam of 38.92 keV (A = 0.3184 A). Pul- 
verized sample was filled into a 1 mm diameter cylindrical 
Kapton capillary and the data collection was carried out 
in a forward scattering geometry using a Perkin Elmer 2- 
D detector mounted orthogonal to the beam path 378.3 
mm away from the sample. 

Single crystals of /3-FeSe were also mounted on glass 
fibers for examination using an Oxford - Diffraction 
Xcalibur 2 CCD 4-circle diffractometer with graphite 
monochromated MoKa radiation. Elemental and mi- 
crostructure analysis were performed on several /3-FeSe 



crystals as well as on the particular crystals chosen for 
resistivity and magnetization using energy-dispersive X- 
ray spectroscopy in a JEOL JSM-6500 scanning electron 
microscope (SEM). 

Sample dimensions were measured with an optical mi- 
croscope Nikon SMZ-800 with 10 /im resolution and M/H 
values were corrected for straw background at each (T,H) 
of the measurement, real sample volume and demagneti- 
zation factor. Thin Pt wires were attached to electrical 
contacts made of Epotek H20E silver epoxy for a stan- 
dard four-probe measurement with current flowing in the 
(101) plane of the tetragonal structure. Magnetization 
and resistivity measurements were carried out in a Quan- 
tum Design MPMS and PPMS respectively. 



III. RESULTS 

Fig. 1 shows Fe-Se and LiCl-CsCl phase diagrams J^ii^ 
The superconducting PbO-type /3-FeSe is a low tempera- 
ture crystallographic phase that decomposes into Fe and 
hexagonal NiAs phase (a-FeSe) at 457° C (Fig. 1(a)). It 
coexists with hexagonal a-FerSeg below 300° C for cer- 
tain Fe-Se stoichiometry.^^ Previous attempts (for exam- 
ple Refs. 10, 11, 16) to prepare single crystals of /3-FeSe 
involved nucleation and growth using KCl/NaCl flux or 
vapor transport reactions. We choose a LiCl - CsCl flux 
method of synthesis due to the presence of a low temper- 
ature eutectic at 326° C, well below the decomposition 
temperature of /3-FeSe.^^ The crystal growth possibly in- 
cludes nucleation of FerSes above 700° C and structural 
phase transition at low temperatures. As opposed to 
crystals grown in KCl,— the low temperature eutectic 
(Fig. 1(b)) allows for complete transition to tetragonal 
/3-FeSe from (457 - 300) °C in a large fraction of crystals 
grown in a batch. Platelike FeSe crystals with the (101) 
plane exposed and elongated in one direction up to 1.5 x 
0.5 X 0.05 mm^ can be separated by dissolving the flux 
in de-ionized water and rinsing in ethanol. 

Crystals of /3-FeSe were separated from a-FerSeg by a 
combination of a permanent magnet and a powder X - ray 
diffraction (XRD) spectra from a Rigaku Miniflex with 
CuKa radiation (A = 1.5418 A). We observed only (hOl) 
peaks (Fig. 2) of the tetragonal phase in selected crystals 
for further analysis (Fig. 2 inset). Crystals contaminated 
with FerSeg and/or oxides showed additional peaks in the 
spectra. 

Fig. 3 and Fig. 4 show electron density maps and el- 
emental analysis of as-grown /3-FeSe single crystals. The 
relative stoichiometry of multiple points on the as-grown 
crystals were measured and a composition of Feo.99(4)Se 
was obtained. Electron density maps of these crystals 
confirmed a uniform distribution of Fe and Se. 

Crystals are moderately air sensitive. After one week 
of air exposure, a selenium oxide layer is visible on 
the crystal surface. After about a month of air expo- 
sure Fe304 is detectable in laboratory X-ray diffraction 
(XRD) and by the observation of a Verwey transition 
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FIG. 1. (Color online) Fe-Se (a) and CsCl - LiCl phase dia- 
grams (b). The presence of the low temperature eutectics (b) 
enables long annealing in the liquid below 457° C and com- 
plete transition from a-FeSe to /3-FeSe in the large fraction of 
the crystals in a batch. 



in M(T).^^ The first attempts to carry out synchrotron 
XRD experiments on samples that were exposed to air 
for several months revealed the presence of multiple ad- 
ditional phases including appreciable amounts of sev- 
eral selenium oxide and iron oxide phases. It has been 
found that while pure stoichiometric FeSe crystals can be 
grown, these tend to degrade through oxidation over the 
course of time. This suggests that the surface of FeSe 
crystals may be Se-terminated and that selenium oxide 
forms first, with further degradation involving iron oxide 
phases or Fe2[Se03]2 as well. Results reported here are 
based on samples whose exposure to air was minimized, 
and only traces of selenium oxide were found in crystals 
that were pulverized for the powder XRD experiments. 

Synchrotron powder XRD data of FeSe sample were 
successfully refined using a two phase structural model 
(Fig. 5). The best fit contained 96.1 mol% (90.7 wt%) 
of FeSe of P4/nmm space group, with a = 3.7622(2) A, 
c = 5.5018(5) A, with Se at (1/4,1/4,0.2624(1)) and Fe 
at (3/4,1/4,0). Compared to the high Tc stoichiomet- 
ric polycrystalline /3-FeSe,^^^^^ the unit cell parameters 
are reduced by 0.3% (a -axis) or 0.4% (c-axis), whereas 
c/a is smaller or identical. The anisotropic atomic dis- 
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FIG. 2. (Color online) Powder X-ray diffraction spectra on a 
/3-FeSe single crystal shows no impurity phases and (hOl) crys- 
tal orientation at 300 K. The data were shown by (+), the fit 
is given by the top solid line and the difference curve (bottom 
solid line) is offset for clarity. Allowed crystallographic re- 
flections are given as vertical tick marks for /3-FeSe (top line) 
and FerSes (bottom line). Bragg peak (201) for 2Q — 51.1 
clearly distinguishes /3-FeSe from FerSes, in addition to mag- 
netic properties. Inset shows typical /3-FeSe single crystals. 
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FIG. 3. (Color online) Microprobe electron density maps of 
as- grown /3-FeSe single crystals. 



placement parameters (ADP) ratio U33/U11 is the ratio 
of thermal vibrations along crystallographic c and a axes 
in tetragonal strucutre. The U33/U11 was 1.12 for Se 
and 1.41 for Fe. To illustrate this, anisotropic ADP's are 
shown as exaggerated thermal ellipsoids in Fig. 5 inset. 
In Van der Waals bonded crystals, such as FeSe, U33/U11 
ADP ratio is expected to be larger than 1 and the ob- 
served ratios are within the expected range. A somewhat 
larger ADP ratio of Fe, 1.41, suggests that it is under- 
bonded and can move along the c-axis. The FeSe4 units 
are found to deviate from perfect tetrahedra, with an 
Fe-Se distance of 2.379(5) A, and tetrahedral angles of 
104.5(5) and 112.0(5) degrees. The anisotropic ADP ra- 
tio (tetrahedral angles) are smaller (equal) than the val- 
ues obtained for /3-FeSe polycrystals on powders contain- 
ing several magnetic impurity phases. ^^^^^^^^However, the 
anisotropic ADP ratio observed here is similar to val- 
ues found in pure Fei.osTe.^^ As expected, the observed 
tetrahedral bond angles deviate from the ideal tetrahe- 
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FIG. 4. (Color online) Elemental analysis of as-grown /3-FeSe 
single crystals. 
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FIG. 5. (Color online) Structural refinement of FeSe syn- 
chrotron powder X-ray diffraction data over a narrow range of 
scattering angle, 2^, taken at 300 K. Background subtracted 
data are shown as (+), fit is given as a top solid line, and 
the difference curve (bottom solid line) is offset for clarity. 
Allowed crystallographic reflections are given as vertical tick 
marks: /3-FeSe phase (bottom), and SesOs impurity phase 
(top) due to sample oxidation. Inset shows /3-FeSe structure. 



dral angle found in iron based superconductors with op- 
timal Tc.^^ The second phase, constituting 3.9 mol% (9.3 
wt%), was found to be SeaOg with Pmc21 space group 
with the refined lattice parameters a = 4.977(1) A, b = 
4.388(2) A, c = 15.377(2) A. No other phases were ob- 
served. Within the main phase we investigated in detail 
the issue of stoichiometry and occupancy of the intersti- 
tial site, Fe(2) at (3/4,3/4,z) (Fig. 5 inset). The stoi- 
chiometry was found to be Fei. 00(2) Sei. 00(3) • 

Difference Fourier analysis (DFA) is a standard 
method to find missing electron densities in refined 
atomic structures. In this technique, the difference be- 
tween the observed and the calculated (model based) 
Fourier maps is used to locate missing atoms in atomic 
structures. In this study, DFA did not reveal any appre- 
ciable electron density at the interstitial Fe(2) positions. 
Attempts to exphcitly refine Fe(2) site occupancy yielded 
0.00(1), in agreement with the DFA and strongly sug- 
gesting that no iron resides on this site. However, DFA 
indicated a possibility for additional electron density in 




50 100 150 200 250 300 
T(K) 

FIG. 6. (Color online) Temperature dependence of the re- 
sistivity for current in (101) plane with Tc onset and zero 
resistance of 12.0(1) and 9.2 (2) K, respectively. Inset shows 
precession pattern of the (101) plane of the same single crys- 
tal. All spots can be indexed within /3-FeSe space group with 
no impurities present. The large mosaic is visible but the 
impurity free unit cell parameters are in agreement with pub- 
lished (see text). 



the vicinity of Se. This, along with the observation of rel- 
atively large anisotropic ADPs (U33) of Se and Fe, may 
point to the presence of static and/or dynamic disorder 
associated with these sites. A small number of Se vacan- 
cies may lead to relaxation of the surrounding Fe atoms, 
resulting in static and/or dynamic disorder. 

The resistivity p(T) of LiCl/CsCl flux grown crystals 
smoothly changes to linear at low temperature where the 
onset of Tc and zero resistivity were observed at tempera- 
tures about 1 K or more higher than in polycrystals (Fig. 
6).^ The residual resistivity ratio (RRR) of 14 indicates 
good crystal quality whereas the single crystal diffraction 
pattern shows no impurities present. Reciprocal space 
planes (hkO) and (hOl) were reconstructed from several 
series of CCD frames (inset in Fig. 6). A mosaic struc- 
ture is observed perpendicular to the c - axis, consistent 
with the arrangement of FeSe layers in the structure. The 
observed X-ray reflections are all consistent with the /3- 
FeSe structure. M/H exhibits weak temperature depen- 
dence for both H ± (101) and H||(101) (Fig. 7(a)). Below 
135 K the M/H signal drops and then remains constant 
below the structural transition temperature ~ 100 K^ 
This is more pronounced for H ± (101). Low temper- 
ature M/H taken in H = 10 Oe confirms superconduc- 
tivity (Fig. 7(b)) below 9.0(2) K. Extrapolation of Airx 
data to T = gives about 60 % of diamagnetic screen- 
ing. Complete p transition and partial superconducting 
volume fraction have been observed in SmFeAsOi-a^F^^ 
and CaFei_a;Coa;AsF,^'^^ where temperature dependent 
magnetic moment coexists and inversely scales with the 
superconducting volume fraction. 

In /3-Fei.oiSe polycrystals contaminated with magnetic 
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FIG. 7. (Color online) Anisotropy in high temperature M/H 
for H = 10 kOe and superconducting volume fraction for H 
± (101) measured in 10 Oe field. 



impurity phases a static moment was found above 1 
GPa,^^ was ascribed to traces of Fe impurity at am- 
bient pressure due to its weak nature^*^ or was not 
detected. ^'^^ In our crystals, s - shape of M(H) for 
H||(101) (Fig. 8(a) is typical of a type-II superconduc- 
tor with a superimposed isotropic weak ferromagnetic 
(WFM) moment both above and below Tc.^^ In contrast, 
the M(H) curves are symmetric for H ±(101) at T = 1.8 
K (Fig. 8(b) within the experimental resolution (0.05 
emu/cm^) with no evidence for WFM below Tc. 



IV. DISCUSSION 

Is the observed WFM intrinsic or extrinsic? Fig. 8(c) 
shows background-subtracted data compared to a simu- 
lated pattern for 1 mole % of commonly found magnetic 
phases. Based on the scattering power and distribution 
of peaks in our synchrotron powder XRD we can exclude 
contamination by a-FeSe and FerSeg (Fig. 8(c) and Ap- 
pendix A), leaving only elemental Fe, or some unknown 
ferromagnetic phase with lattice periodicity commensu- 
rate with /3-FeSe^ as a possible source of magnetic con- 
tamination. By focusing on information sensitive to the 
total number of Fe ions in a given volume we make several 
observations that support the intrinsic WFM scenario in 
as-grown crystals. Mossbauer quadrupole splittings and 
isomer shifts for binary Fe-Se materials differ by up to 
a factor of 4 whereas the signal component originating 
in the /3-FeSe phase dominates the Mossbauer spectrum 
below room temperature in 50 kOe even for the sample 
containing more than 6 mole % of a-FeSe, Fe304 and ele- 
mental Fe impurities. ^^^^^^^^ Moreover, the coercive field 
of bulk Fe or Fe nanoparticles (/ioHc < (0.4 - 2.5) kOe)is 
several orders of magnitude lower than in Fig. 8(a,b) 
(15 kOe)^i^ Furthermore we observe the signature of 
the structural phase transition in the normal state (Fig. 
7(a)), implying that a considerable fraction of the M(T) 
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FIG. 8. (Color online) Magnetization isotherms below and 
above Tc for magnetic field parallel (H||(101)) (a) and per- 
pendicular to (101) plane (H ± (101)) (b). Right hand side 
scale shows magnetic moment assuming molar mass of /3-FeSe 
(c) Comparison of background - subtracted diffraction data 
with simulated impurity phases at 300 K. Selected 20 region 
is shown where simulated impurity peaks are clearly visible. 



signal must come from the /3-FeSe. Finally, WFM due to 
an unknown Fe - based high Tc ferromagnetic phase^^ is 
unlikely to have /ioHc in kOe range and is expected to 
provide an isotropic constant (or increasing) background 
(bias) to M(H) loops both above and below Tc, as in ref. 
11. This is in contrast to anisotropic M(H) below Tc (Fig 
8(a,b)). Note that M(H) in Fig. 8(b) is symmetric with 
respect to the M = line. This suggests that most of 
WFM signal in /3-FeSe crystals is unrelated to extrinsic 
impurities. The WFM signal is well reproduced in sev- 
eral crystals grown in one batch and in crystals grown 
from multiple batches, whereas the magnitude of WFM 
increases in time (Appendix B). 

The proposed spin-density- wave (SDW) mechanism 
of magnetic order may not apply to iron chalcogenides 
and perhaps it could be more complex even for iron 
pnictides.^^^^^ Density functional calculations indicate 
that magnetic state with 0.15 /i^/mole is induced in 
Fei_a,Se for x = 0.0625.^^ The saturation moment \Ms\ 
observed (Fig 8(a,b)) is about 1/150 smaller and would 
correspond to an Fe deficiency of less than 0.4 atomic %, 
a value below resolution limit of our diffraction measure- 
ment. However, in this calculation^^ the partial density 



of states (PDOS) of Fe(2) dominates the total density of 
states (DOS) and, more importantly, the stoichiometry 
variation. Negligible occupancy of Fe(2) within experi- 
mental error in our crystals implies that Fe nonstoichiom- 
etry is not the dominant mechanism of WFM. On the 
other hand, it is possible that some fraction of the WFM 
arises due to a Se vacancy induced magnetic cluster^ 
The Se - Se distances are Van der Waals distances and 
may produce Se - Se time dependent bonding. We cannot 
distinguish between static and dynamic displacements 
but since the refinement results are giving 1:1 stoichiom- 
etry, then vacancies could be equally distributed on both 
sites. In particular, theory predicts that the main effect 
of Se displacement would be to shift Fe(l) towards the 
vacant site, shifting the Fermi level E^? into a sharp peak 
in the DOS that would promote a more stable magnetic 
state than in material without Se defect s.^^ The net mo- 
ment at 1.25 mole % Se deficiency is expected to be in 
10~^ /i^/niole range. This is arising from both Fe(l) and 
Fe(2) contributions. Since the theoretical contribution of 
Fe(2) PDOS at the Fermi level is about 50% of the total 
DOS,^^ the calculated moment is somewhat higher but 
generally in line with the observed \Ms\ ^ (1.0±0.5)-10~^ 
/i^/mole above Tc in our crystals with neglibible Fe(2) 
occupancy. This is different from most FeAs supercon- 
ductors where small moment magnetic order from a SDW 
mechanism is found below the structural transition. This 
is also different from Fei+^yTe where subtle crystal chem- 
ical effects, with both Fe(l) and Fe(2) occupied, induces 
WFM and structural and magnetic differences below the 
magnetostructural transition at 75 K - 55 Ki Since lat- 
tice distortions were also recently found^^ to induce both 
superconducting and magnetic phases in SrFe2As2, this 
suggests that nanoscale defects and short range struc- 
tural features are important in a wider class of iron based 
superconductors. Indeed, there is emerging evidence that 
both conducting and magnetic properties in the recently 
discovered KxFe2-ySe2 superconductors are governed by 
Fe vacancies ^i^ In /3-FeSe, defect - induced magnetism 
coexists with superconductivity that sets in far below 
Tc. Though rather unlikely, our analysis allows for some 
contribution of different impurity phase to WFM. The 
unknown high temperature FM phase would be present 
in quantities too small to detect by diffraction and/or 
would have the lattice periodicity commensurate with 
/3-FeSe and its moment would anisotropically diminish 
between 30 K and 1.8 K. 



In summary, single phase superconducting single crys- 
tals of /3-FeSe have been synthesized. Unlike isostruc- 
tural Fei+^Te, the Fe(2) site is not occupied at all in 
these samples. The ADP anisotropy is consistent with 
dynamic disorder/defects associated with Fe and Se sites 
and/or Se vacancies. We present evidence for intrinsic de- 
fect induced WFM which anisotropically diminishes with 
an increase in the superconducting volume fraction. 
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VI. APPENDIX A: MAGNETIC IMPURITIES 
FROM DIFFRACTION 

Magnetization measurements were performed on a 
sample volume V = 4.532-10"^ cm^ (a = 0.223476 cm, 
b = 0.0555613 cm, c = 0.00365025 cm), corresponding 
to (1.45±0.23)-10"^ emu MPMS signal. Assuming that 
magnetization (M) of /3-FeSe is negligible when compared 
to impurity magnetization, we discuss possible magnetic 
impurity levels. 

First we assume that the sample contains 1 mole % of 
elemental a-Fe since any higher Fe content would have 
been detected (Fig. 8(c)). In order to obtain volume ratio 
we need to multiply mole (i.e. formula unit) ratios of Fe 
and /3-FeSe by (M/D) where M is the molar mass and D is 
the density in g/cm^. Using M(Fe) = 55.8 g/mole, D(Fe) 
= 7.83 g/cm^ M(/3-FeSe) = 134.8 g/mole, D(/3-FeSe) = 
5.72 g/mole, we get 0.3 volume %. 



Imole(Fe) • 7.12^^ 

^ ^ mole 

100mole{/3 - FeSe) • 23.56^ 



= 0.003 



(1) 



That would correspond to 1.37-10" cm'^ volume of 
Fe in our sample. What would be its magnetic contri- 
bution? By dividing D/M we get 0.14 mole/cm^ of Fe. 
In 1.37-10"^ cm^ iron volume we have 1.92-10"^ mole 
Fe. Iron saturation moment is 2.2 /iB/mole, therefore 
that volume would have M = 4.22-10"^ jib- Using con- 
version factor 5585 emu//iB we obtain M = 2.36-10"^ 
emu. This is more than 100 % of raw M (emu) signal in 
MPMS. Hence, it is possible that Fe impurity content in 
our crystal (mole % and consequently volume %) is below 
detection capacity of synchrotron powder X-ray diffrac- 
tion data since only a fraction of 1 mole % Fe would 
generate such signal. This is about 0.05 atomic % of Fe. 

Similarly, 1 mole % impurity of a-FeSe is approxi- 
mately identical to 1 % volume since a-FeSe has the iden- 
tical molar mass M and 95% of /3-FeSe density. Therefore 
1 % mole a-FeSe impurity would correspond to 4.5-10"^ 
cm^ volume (using eq. (1)). Applying the same argument 
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FIG. 9. (Color online) Structural refinement of FeSe syn- 
chrotron powder X-ray diffraction data over the full 20 range 
taken at 300 K. Data (background subtracted) are shown 
as (+), fit is given as a top solid line, and the difference 
curve (bottom solid line) is offset for clarity. Bottom verti- 
cal tick marks represent reflections in the main /3-FeSe phase 
{P4/nmm) , while top tick marks denote reflections in SesOs 
{Pmc21). 



as above, 4.5-10"^ cm^ volume a-FeSe contains 1.9-10"^ 
mole a- FeSe. The a- FeSe saturation moment is Mg '^0.2 
/iB/mole4^ Therefore such volume would contribute M 
= 3.81-10"^ /J.B' Using conversion factor 5585 emu/ /j^b 
we obtain M = 2.13-10"^ emu. Our raw MPMS signal is 
70 % that value. However, 0.7 mole % of a- FeSe would 
have been detected, if present (70 % of its peak height, 
(Fig. 8(c), main text)). 

Finally, 1 mole % of FeySeg impurity (using M(Fe7Se8) 
= 1022.59 g/mole, D(Fe7Se8) = 6.43 g/cm^ and eq. (1)) 
would correspond to 6.75 % of measured sample volume, 
which is 3.059-10"^ cm^. Since the expected FerSeg sat- 
uration moment is M^ ~ 80 (emu/cm^ )^^ we would ex- 
pect that such volume would contribute with 2.44- 10~^ 
/j^B- Our raw MPMS emu signal is 6 % of that value, but 
still above the threshold of scattering power detection in 
synchrotron experiment. If we multiply the observed in- 
tensity of 2.92 A~^ FerSes peak (Fig. 8(c) main text) by 
0.06, it is still above the background. 



VII. APPENDIX B: REPRODUCIBILITY, 

IMPURITIES AND MAGNETIC SIGNAL OVER 

TIME 

Figure 9 shows Rietveld refinement over the full 20 
range with only selenium oxide present due to oxidation. 

Powder X-ray diffraction taken on single crystals found 
no extrinsic peaks in 10 out of 10 crystals that were sep- 
arated by a magnet. Only clean /3-FeSe single crystals 
are used in further analysis. Samples that were contam- 
inated by FeySeg (from synthesis or oxidation) showed 
dominant FM hysteresis loop below Tc with only traces 
of type-II superconductor MvsH. Resistivity and magne- 
tization data were well reproduced in several indepen- 



0.0 



-0.2 



-0.4 






-0.6 



-0.8 




0123456789 10 
T(K) 

FIG. 10. (Color online) An example of superconducting Tc 
of independently grown crystal taken in 10 Oe field. 




Hl(lOl) plane 

sample 2 
— o— 1.8K 
— □— 4K 
— A— 30K 
— V— 250K 




-40 -20 20 40 
H (kOe) 



-40 -20 20 40 
H (kOe) 



FIG. 11. (Color online) (a) Sample 2 M/H above and below 
Te for H||(101) near M = 0. (b) Sample 2 M/H above and 
below Tc for H±(101) near M = 0. 



dently grown crystals from multiple batches. Both bulk 
Tc, as measured by Airx and resistive Tc (as defined in 



the text) varied by ATc = ± 1 K. This variation is prob- 
ably due to sample degradation induced by variable air 
exposure. However, the onset of resistive Tc was always 
above the bulk. Fig. 10 shows example of superconduct- 
ing Tc and volume fraction for different, independently 
grown sample. Figures 11(a) and 11(b) show hysteresis 
curves below Tc for independently grown crystal from the 
same batch (sample 2) as the crystal used in the main 
text (sample 1). Unlike sample 1 (that was measured 
within a day from the moment of its synthesis), sample 
2 was exposed to air for several days. In addition, it was 
slightly heated when sealing in quartz tube and kept in 
the low vacuum (several Torr) for about 3 months. Dom- 
inant type-II superconducting MvsH hysteresis is evident 
whereas small WFM is superimposed on the main signal. 
As expected, the magnitude of WFM signal is about two 
times larger in sample 2 than in crystal that was exposed 
to air for shorter time (Fig. 8(a,b), main text). Based on 
our synchrotron powder X-ray results taken on samples 
that were exposed to air for several months, the larger 
WFM magnitude should also originate from Fe - based 
compounds (impurities) that form over a course of time 
(see main text). The thickness of sample 2 was identi- 
cal and ab plane was considerably smaller than sample 
1, hence the crystal was more isotropic (a = 0.1535 cm, 
b = 0.036 cm, c = 0.0036 cm). Larger irreversibility 
fields are expected when geometric edge barrier for vor- 
tex penetration dwarfs pinning at the inhomogeneities of 
the material. If intrinsic, larger WFM signal in sample 
2 is expected to originate from more defects that would 
cause larger irreversibility field if the bulk pinning on 
inhomogeneities is dominant. However, for thin super- 
conducting strips geometric (specimen shape dependent) 
barrier is dominant and larger irreversibility fields are 
expected for more anisotropic samples.^^ Both sample 1 
and sample 2 are rather thin (c/a = 0.016, c/b = 0.06 
for sample 1 and c/a = 0.023, c/b = 0.1 for sample 2), 
whereas sample 1 is more anisotropic and is expected to 
have larger irreversibility field. 
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